arXiv:1505.06884vl [gr-qc] 26 May 2015 


The flyby anomaly: 

A case for strong gravitomagnetism ? 

L. Acedo* 

Institute Universitario de Matematica Multidisciplinar 
Building 8G, 2° Floor, Camino de Vera, 
Universitat Politecnica de Valencia, 46022, Valencia, Spain 
e-mail: luiacrod@imm.upv.es 

May 27, 2015 


Abstract 

In the last two decades an anomalous variation in the asymptotic 
velocity of spacecraft performing a flyby manoeuvre around Earth has 
been discovered through careful Doppler tracking and orbital analysis. 

No viable hypothesis for a conventional explanation of this effect has 
been proposed and its origin remains unexplained. In this paper we 
discuss a strong transversal component of the gravitomagnetic field as 
a possible source of the flyby anomaly. We show that the perturba¬ 
tions induced by such a field could fit the anomalies both in sign and 
order of magnitude. But, although the secular contributions to the 
Gravity Probe B experimental results and the Lense-Thirring effect in 
geodynamics satellites can be made null, the detailed orbital evolution 
is easily in conflict with such an enhanced gravitomagnetic effect. 

Keywords: Flyby anomaly. Extended General Relativity, Celestial me¬ 
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1 Introduction 


In the last decades we have entered an era of high-precision astrodynam- 
ics in which spacecraft, planets and satellites can be tracked with radar 
ranging, laser ranging and Doppler techniques [Downs and Reichley(1975), 
Williams and Dickey(2002), Duev et al.(2012)]. These improved methods 
have allowed the testing of orbital trajectories in General Relativity to un¬ 
precedented accuracy. In particular, the determination of the (3 and 7 pa¬ 
rameters in parametrized post-Newtonian dynamics is now possible to the 
level of ~ 10“®, ruling out alternative models, such as Brans-Dicke theory, 
or putting stringent constraints upon them [Perivolaropoulos(2010)]. 

Error checking and detailed modelling of spacecraft geometry and thermal 
properties has also proven useful in the improvement of accurate navigation 
methods. A finite-element model for the anisotropic thermal radiation from 
the on-board radioisotope generators has been build for the Pioneer space¬ 
craft. This model has explained away the anomalous recoil acceleration di¬ 
rected toward the Sun which has been found both for Pioneer 10 and Pioneer 
11 spacecraft [Rievers and Lannnerzahl(2011), Tnryshev et al.(2012)]. 

Another anomalous behaviour of spacecraft was reported in 2008 by An¬ 
derson et al [Anderson et al.(2008)]. These authors have analyzed the data 
for several Earth flybys of deep-space missions that took place between De¬ 
cember 1990 and September 2005. Flybys are a common manoeuvre in space¬ 
craft missions which allows the spacecraft to gain or lose of heliocentric energy 
with the purpose of reaching their objective [Borcherds and McCauley(1994)]. 
An analysis of the data for these flybys have shown X-band Doppler residuals 
that are interpreted in terms of a change of the hyperbolic excess velocity, 
Poo) of a few mm/s. Anderson et al. have proposed the phenomenological 
formula: 

= K{cos 6 i - cos 60 ) , (1) 

^00 

where 6 i, 60 are the declinations for the incoming and outgoing osculating 
velocity vectors and iP is a constant. The value of K seems to be close to 
2ueRe/c, where cug is the angular rotational velocity of the Earth, is the 
Earth radius and c is the speed of light. Although this formula works reason¬ 
ably well for the six flybys studied in the paper, the proposal for the relation 
of K with the Earth’s tangential velocity at the Equator is an unjustified and 
speculative hypothesis, taking into account that the flybys of other planets 
with different rotational velocities and radii have not been considered. 
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In the aforementioned paper, Rievers et al. [Rievers and Lanimerzahl(2011)] 
also analyzed the thermal radiation pressnre on a model of the Rosetta space¬ 
craft and found it insufficiently small to account for the anomaly as well as 
wrong in sign. Other conventional effects that have been considered and es¬ 
timated are: atmospheric drag, ocean or solid Earth tides, charge and mag¬ 
netic moment of the spacecraft. Earth albedo. Solar wind and spin-rotation 
coupling [Lannnerzahl et ah (2008)]. lorio has also computed the effect of 
the perturbation induced by the Earth’s gravitomagnetic held as predicted 
by General Relativity [Iorio(2009)]. This effect is hve orders of magnitude 
smaller than the observed velocity change. We will see that the enhanced 
gravitomagnetic effect, proposed from a phenomenological point of view in 
this paper, is modulated by a parameter f3 much larger than the prefactor 
rs/r (the quotient of the Schwarzschild radius and the distance from the 
planet’s center to the spacecraft) of the standard gravitomagnetism of a ro¬ 
tating planet in General Relativity. Another attempt to explain the anomaly 
uses the transversal Doppler effect in Special Relativity [Mbelek(2009)] but 
this explanation is also dismissed because the spacecraft movement is not 
always transversal to the Earth surface and the anomaly appears also in the 
ranging data as well as in the Doppler data. It has also been suggested that 
a second-order error in the integration method could manifest as the alleged 
anomaly [Acedo(2014)]. 

Moreover, a mismodelling of the coordinate systems may also be involved 
in the anomaly because the fitting equation of Anderson et al. [Anderson et al.(2008)] 
involves only geometrical parameters [Lanimerzahl(2011)]. The subtleties of 
the energy transfer process during planetary flybys have also been analyzed 
in detail for some real missions in the hope of gaining some insight on the 
origin of the flyby anomaly but to no avail [Anderson et al.(2007)]. In the 
flyby manoeuvres the spacecraft gains or lose energy with respect to the 
barycenter of the Solar System but the strange fact about the anomaly is 
the apparent energy change observed in the Earth-centered system. In this 
paper we will study transient variations of the perturbed velocities and posi¬ 
tions for the spacecraft as the consequence of a two-body interaction among 
the spacecraft and a phenomenological proposal for a gravitomagnetic field 
of the Earth. So, our results are expected not to be affected by the definition 
of the coordinate frame. 

Some proposals for a non-conventional explanation have also been made: 

Adler has worked out a model for an halo of dark matter particles bound by 
the Earth gravitational held. The explanation of the anomaly requires very 
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strong constraints on the cross sections, the energy of the particles and the 
dominance of inelastic and elastic processes along different trajectories to ac¬ 
count for the increase or decrease of the spacecraft velocities [Adler(2009)]. 
A hve dimensional theory of gravity has also been considered by Gerrard 
and Sumner [Gerrard and Sumner(2010)]. Other proposals include radical 
modihcations of well established principles such as the principle of equiv¬ 
alence [McGulloch(2008)] or Lorenz invariance [Gahill(2008)]. Busack has 
also derived a phenomenological formula for a velocity dependent potential 
with three adjustable parameters which hts the observed discrepancies in 
the flybys [Busack(2013)]. However, these radical proposals are unlikely to 
receive further support because they are disconnected from the theoretical 
paradigms or other domains of empirical data. 

In this paper we take the data for the flyby anomalies at face value and we 
consider a scenario in which the anomalies are the consequence of a modelling 
force unpredicted by the current theory of gravitation. To select this model 
we follow an intuition given by the authors of the seminal paper on this topic 
[Anderson et al.(2008)]. In their work, Anderson et ah suggest that the 
latitude dependence of the anomaly could be related with a frame-dragging 
effect much larger than the one predicted by standard General Relativity 
[lorio et ah (2011)]. It is well-known that linearized General Relativity can 
be interpreted in terms of a gravitoelectric held (containing Newton’s law) 
and a gravitomagnetic held [Rindler(2006)], which for a rotating planet takes 
the form: 

B = -i^((^)'in-3(n.f)f] . (2) 

where rs = 2GM/c^ is the Schwarzschild radius of the body, R is the plan¬ 
etary radius and hi is the angular velocity vector. This held imparts an 
acceleration in a test particle given by a Lorentz force law: 

2 

a = —V X B , (3) 

c 

which justihes the denomination of B as gravitomagnetic held. This ehect, 
jointly with direct detection of gravitational waves, is one the most feeble 
and difhcult predictions of General Relativity to be measured but, after sev¬ 
eral decades of planning, the gravitomagnetic held of Earth (as well as the 
larger geodetic ehect due to the movement in the curved space-time around 
the planet) has hnally been studied experimentally and the hnal result is 
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in good agreement with the theory [Everitt et ah (2011)]. The frame drag¬ 
ging drift was measured with four high precision electrically suspended gy¬ 
roscopes mounted in the Gravity Probe B satellite in a polar circular orbit 
642 km height. Being a secular effect, it piled up orbit after orbit during 
one year of data collection and a hnal effect of —3.72(0.72) x 10“^ arcsec 
was derived. The measured effect is in good agreement with the predic¬ 
tions of General Relativity but the 20% error is still high. It is remarkable 
that this data was obtained after a careful modelling of a much larger effect 
on the alignment of the gyroscopes caused by electrostatic interactions be¬ 
tween the rotors and their housing [Everitt et al.(2009), Reiser et al.(2009), 
Muhlfelder et ah (2009), Silbergleit et ah (2009)]. 

The prefactor r^/r in Eq. (2) is ~ 1.26 x 10“® for the Gravity Probe B 
orbit and this gives us an idea of the small magnitude of the frame-dragging 
effect. One could think that the results of this mission rule out the possibility 
of any other gravitomagnetic effect beyond conventional General Relativity. 
However, we must take into account that Gravity Probe B was not designed 
to detect non-secular contributions. A contribution of such a kind could 
arise from a gravitomagnetic held tangential, at every point, to the celestial 
parallel in the form: 

B = R(r,0)0, (4) 

where 6 is the polar angle and 0 is the azimuthal angle. Notice that for any 
circular orbit around the Earth there is always a symmetric point, for any 
other point, in which the vector held in Eq. (4) has opposite direction to the 
vector in the hrst point of the trajectory. This implies that the gyroscope 
precession induced by this held do not add up orbit after orbit but cancels out. 
Moreover, we will calculate the contributions to the Lense-Thirring precession 
of the longitude of the ascending node of the geodynamics satellites to check 
the viability of our proposal in the context of the most accurate measurements 
of this ehect to date[Iorio et al.(2011)]. 

On the other hand, the ehects of this new held would manifest on highly 
eccentric elliptical orbits or in asymmetrical hybys. This could explain why 
the anomaly discovered in the hybys of the Earth only appears on these 
particular orbits. Our objective is to propose a reasonable form for B{r, 9) in 
Eq. (4) by dimensional and continuity arguments and analyze its potential 
as an explanation of the hyby anomaly. 

The paper is organized as follows: In section 2 we discuss the form of 
the transversal gravitomagnetic held and develop a perturbative approach 
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to analyze its influence on flybys. Section 3 is devoted to the calculations 
of the perturbations for the NEAR, Galileo, Cassini, Rosetta flybys. We 
also give a prediction for the most recent Juno flyby. The effects on the 
Gravity Probe B experiment and the LAGEOS, LARES and GRACE Lense- 
Thirring precession is discussed in Section 4. Finally, we remark on the 
context and further development of the idea of a transversal gravitomagnetic 
held in Section 5. 


2 Flybys in a transversal gravitomagnetic field 


We should begin by proposing a reasonable hypothesis for a transversal com¬ 
ponent of a gravitomagnetic held. Firstly, we can expect it to be proportional 
to the total angular momentum of the Earth. Moreover, the continuity con¬ 
dition B{r,6 = 0) = B{r, 6 = it) =0 must be verihed. We can also argue 
that B{r, 6) is an even function of the polar angle, B{r, —9) = B{r, 6). From 
these conditions alone we suggest that a simple form for a transversal grav¬ 
itomagnetic held is: 

B(r, 6^) =sin6* cos6* 0 , (5) 

r 

where is the angular velocity of the Earth, R is its radius and /S is a 
constant to be determined. The reason for the particular dependence in the 
polar angle, 9, and the decrease as the inverse of the distance to the center 
of the Earth, r, will be revealed in Sec. 4 in relation with the constraints 
imposed by other precise measurements of the Earth’s gravitational held. 
The unit azimuthal vector, cf) can easily be expressed in the cartesian celestial 
coordinate system whose x axis points toward the hrst point of Aries (the 
point where the Sun crosses the celestial equator during the Vernal equinox). 
In this system the axis points towards the celestial pole and this dehnes 
the remaining y axis in a right handed coordinate system. So, we have: 


sin 9 cos (j) 
sin 9 sin (j) 


r ■ i 

|r| 

r j 

|r| 


( 6 ) 


In order to calculate the contribution of the small acceleration a = v x B to 
the perturbation of the hyperbolic orbit in a hyby we must take into account 
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also the additional second-order effects induced by the Earth’s Newtonian 
gravitational field, the Sun and the Moon. We must start with some standard 
relations for the celestial mechanics of a flyby hyperbolic orbit. 


2.1 Ideal flyby orbits 


In Fig. 1 we have plotted the trajectory of the NEAR spacecraft during its 
flyby of January 23, 1998. Parameters for this and other flyby orbits are 
given in Table 1. It is expedient to parametrize these hyperbolic trajectories 
in a coordinate system relative to the orbit. So, we choose a unit vector 
along the periapsis direction corresponding to the point of closest approach, 
s, a second unit vector pointing along the direction of the inclination vector 
of the orbit, w, and a third one perpendicular to those two, fi. This third 
unit vector is dehned in such a way that the scalar product with the initial 
radiovector of the spacecraft, Tin, is positive. These vectors are given as 
follows: 

s = cos 9pk + sin 9p cos api + sin 9p sin ap j 
w = cos / k J-sin / cos i J-sin / sin a/j (7) 

h = ±w X s , 


where 9p, ap are the celestial polar angle and right ascension of the periapsis, 
I and are the inclination and the right ascension of the inclination vector 
and the sign in the last expression for h depends on the orientation of the 
orbit. The orthogonal system i, j, k is, obviously, the celestial coordinate 
system. 

In the NEAR flyby case the parameters were (all angles in degrees): 9p = 
57, ap = 280.43, / = 108.0, aj = ap + arccos(—cot/cot6'p) = 358.24. The 
incoming direction is given by 9i = 69.24 and ai = 81.17. In this particular 
case, it can be shown that h = w x s. 

It is also usual to use the eccentric anomaly, r], instead of time, t, or the 
true anomaly, u. The true anomaly is the angle formed by the radiovector 
at any instant and the radiovector at the instant of closest approach to the 
Earth. The relation among t] and u is given by: 

e -f cos u 

coshr^ =-, 

1 J- e cos u 

where e > 1 the eccentricity of the hyperbolic orbit. The time of flight can 
also be given in terms of the eccentric anomaly by: 


t = T(esinh ?7 — rj) , 


(9) 
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Figure 1: Plot of the NEAR flyby orbit (January 23, 1998). The solid vector 
points towards the Sun and the dashed vector points towards the Moon at 
the instant of the closest approach. The asymptotic ingoing direction is also 
shown. 














where the time-scale T = ^ (—a)^//i, a being the semi-major axis of the 
orbit and /i is the product of the gravitational constant and the mass of the 
Earth, /te = 398600.4 km^/s^. 

The equations for the radiovector and the velocity of the spacecraft in the 
ideal hyperbolic orbit are then given by 


r(h) 

v(r7) 


i(cosh 7] — e)s — ay /— 1 sinh rj h 
(sinh T] s 


T(e coshr; — 1) 

— 1 cosh 7] h) , 


( 10 ) 


( 11 ) 


Concerning the tracking of the spacecraft and the inference of the flyby orbit 
from the tracking data there are many practical problems still not analyzed 
in full detail. An accurate flyby model should take into account the most 
recent GGM03 model for the Earth gravitational field as elaborated from the 
data of the GRACE satellite [Tapley et al. (2007)]. This model has achieved 
a large geographical resolution and a precision on the variations of the grav¬ 
itational held at the Earth’s surface down to 1 milligal (10“^ m s“^). On 
the other hand, the position and velocity of the spacecraft is tracked by 
different ground stations around the globe in different relative movements 
with respect to the spacecraft. This should be also be taken into account 
in the modelling of the Doppler effect. Preliminary assessments showed that 
these effects are relatively small and unimportant in connection with the 
anomaly [Lamnierzahl et al.(2008). Paramos and Hechenblaikner(2012)] but 
their rigorous analysis should help in the improvement of the error bars on 
the anomalous accelerations. 


2.2 First and second-order perturbations 

In the following it is useful to use the non-dimensional parameters, r = t/T 
for time and v* = T5w/\a\ for the velocity. A hrst-order perturbation comes 
from the acceleration imparted by the gravitomagnetic field in Eq. (5) as 
follows: 

aMG = 5*(r,6')v* X 0 , (12) 

with B*{r,6) = B{r,6) * T. The cross product in Eq. (12) is readily 
calculated in cartesian coordinates from Eqs. (6), (7) and (11). The position 
of the spacecraft deviates from the ideal hyperbolic orbit discussed in Section 
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2.1 as a consequence of this first-order acceleration. However, this implies 
another second-order contribution from the change of the Earth’s Newtonian 
gravitational field at the new orbital points: 

<^^Earth = 3 /^E ^ ) (^3) 

where /i^ = T/iE/|ap is the scaled gravitational constant times the mass of 
the Earth. As it is well-known, most of the orbital perturbations come from 
the tidal forces exerted by the Sun (and in minor proportion by the Moon). 
If the Sun is located at radiovector R, with the center of the Earth as origin, 
and the spacecraft is found at r the tidal force is given by: 


3-tidal — /^S 



R 


(r2 + R2 _ 2r . R 


, 3/2 


(14) 


with /is = 1.3271244x lO^^km^/s^ as the product of the gravitational constant 
and the mass of the Sun. The change in the position of the spacecraft from 
the transversal gravitomagnetic force induces a second-order contribution 
from the tidal force as follows: 


<^^SUN — /^S" 


(5r 


3/is- 


+ R2 _ 2r • R 
(R — r) • hr 


, 3/2 


(15) 


-|- R2 _ 2r ■ R 


, 5/2 


(R 


Assembling these results in Eqs. (12),(13) and (15) we find that the evolution 
of the perturbations in the position and the velocity of the spacecraft coming 
from the effect of the transversal gravitomagnetic field should be calculated 
from the system: 

hv* 

^MG T ^EARTH T ^SUN • 

By solving this system for the differential equations of motion we should be 
able to compare with the anomalous perturbations of the Earth flybys in the 
next section. To perform this integration it is usually more convenient to 
employ the eccentric anomaly instead of time as defined in Eq. (9). 



d6r* 


dr 

dSv* 


dr 
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Table 1: Parameters for the spacecraft flybys of the Earth. 


Spacecraft 

Date 

e 

a (km) 

^IN 

^OUT 

e. 

/ 

®IN 

0 

NEAR 

1/23/1998 

1.8135 

-8494.87 

69.24° 

161.96° 

57° 

108° 

81.17° 

280 

Galileo I 

12/8/1990 

2.4729 

-4977.24 

77.48° 

124.25° 

64.8° 

142.9° 

86.60° 

319 

Galileo II 

12/8/1992 

2.3194 

-5058.31 

55.74° 

94.87° 

123.8° 

138.7° 

39.47° 

302 

Cassini 

8/18/1999 

5.8525 

-1555.09 

102.92° 

94.99° 

113.5° 

25.4° 

154.33° 

245 

Rosetta 

3/4/2005 

1.3118 

-26710.9 

92.81° 

124.29° 

69.8° 

144.9° 

166.68° 

22. 

Rosetta II 

11/13/2007 

1.5401 

-33417.5 

79.32° 

1-^ 

o 

o 

154.7° 

1-^ 

o 

o 

45.95° 

30^ 

Rosetta III 

11/13/2009 

1.5976 

-25491.1 

108.4° 

65.65° 

97.44° 

155.6° 

31.78° 

27( 

Juno 

9/10/2013 

4.6489 

-3645.92 

104.21° 

50.59° 

123.39° 

47.13° 

215.40° 

344 


3 Numerical results 

In Table 1 we show the parameters for the orbits of several spacecraft flybys 
in which the anomalous variation of the velocity has been found. We have 
displayed only the parameters necessary for the calculations including the 
orbital eccentricity, e, the semi-major axis of the hyperbolic orbit, a, the 
polar and azimuthal angles locating the incoming and outgoing directions, 
the perigee and the inclination vector. These data were collected or inferred 
from the JPL Horizons on-line system data [JPL Horizons(2014)] as well as 
Anderson et ah paper [Anderson et ah (2008)]. We should also know the 
average position of the Sun during the flyby as this is the major source of 
second-order perturbations. These data are tabulated in Table 2 for the eight 
flybys considered. 

Using Eqs. (12), (13) and (15)we have performed a numerical integration 
of the system in Eq. (16) for rj = 0 onwards and backwards. A typical result 
for the Galileo II flyby is plotted in Fig. 2. In this particular case we see that 
the perturbation corresponds to a decrease of the orbital velocity after the 
closest approach to the Earth. The observed variation in the velocity can be 
compared with the differences in the temporal averages of 6V after and before 
this perigee. A simple rough measure of the variation 6V can be obtained as 
the difference among the peak perturbations occurring in the post-encounter 
and pre-encounter regions. In Fig. (3) we have plotted the results for the vari¬ 
ation of 6V in the cases of the NEAR, Rosetta and Galileo II flybys in terms 
of the value of the parameter f3 in Eq. (4). In the NEAR case the largest ve¬ 


il 
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Figure 2: Velocity perturbations induced by the transversal gravitomagnetic 
field in the Galileo II fiyby. In this case we have taken /? = 2.0 x 10“^. Notice 
that the average difference between the sections of the orbits after and before 
the perigee is negative in agreement with the observed decrease during this 
flyby. 
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Table 2: Average distance and average unit vector towards the Sun during 
the flybys in the celestial cartesian coordinate system. 


Spacecraft 

(Rs) (10® km) 

(fs) 

NEAR 

1.4727 

(0.5413,-0.7700,-0.3338) 

Galileo I 

1.4739 

(-0.2594,-0.8852,-0.3838) 

Galileo H 

1.4737 

(-0.2510,-0.8872,-0.3847) 

Gassini 

1.5147 

(-0.8072,0.5403,0.2343) 

Rosetta 

1.4835 

(0.9587,-0.2580,-0.1119) 

Rosetta H 

1.4809 

(-0.6513,-0.6951,-0.3013) 

Rosetta HI 

1.4808 

(-0.6447,-0.7002,-0.3035) 

Juno 

1.4882 

(-0.9615,-0.2479,-0.1075) 


locity increase, amounting to 13.46 mm/s as it was deduced from the X-band 
Doppler residuals, was detected [Anderson et al.(2008)]. This amounted only 
to 1.8 mm/s for the Rosetta flyby. The Galileo II spacecraft approached suf¬ 
ficiently close to the Earth surface (a 303 km altitude at the perigee) for the 
atmospheric drag effect to be important. Nevertheless, Anderson et ah de¬ 
duced an unexpected decrease of —4.6 mm/s. These results are compatible 
with a (3 parameter value in the range 1.4 x 10“^ < /d < 3.0 x 10“^. Similarly, 
in Fig. (4) we plot the results for the Galileo I, Juno and Gassini flybys. For 
this Gassini flyby (whose perigee was much higher for any atmospheric drag 
to be relevant) a decrease of —2 mm/s was derived from the orbital data 
[Anderson et al.(2008)]. Our model can also predict the result for the still 
not analyzed Juno flyby. The prediction is a behaviour very similar to the 
Gassini flyby and, consequently, a similar decrease is expected. On the con¬ 
trary, Anderson et al.’s fitting equation [Anderson et al.(2008), Iorio(2013)] 
predicts an increase around 6 mm/s. However, we must consider that the 
uncertainties could be large in these isolated measurements. 

Finally, we must discuss the Rosetta II and Rosetta III flybys, which have 
been the subject of some polemic [Paramos and Hechenblaikner(2012)]. The 
phenomenological formula by Anderson et ah in Eq. (1) predicts small, but 
measurable, increases of the asymptotic velocity around 0.36 mm/s and 0.46 
mm/s for the Rosetta II and III, respectively. However, a careful analysis 
of the data showed no anomalous behaviour of the fitted trajectory. Gon- 
sequently, despite the success in fitting the previous flyby anomalies, we 
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Figure 3: Modelled anomalous increases of the velocity during the NEAR 
(circles), Rosetta (squares) and the Galileo II flybys (triangles) versus the 
parameter 13. Horizontal solid, dashed and dotted lines correspond to the 
observed values for the NEAR, Rosetta and Galileo II flybys, respectively. 
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Figure 4: The same as Fig. 3 but for the Galileo I (open circles), Juno 
(closed circles) and Cassini (squares) flybys. In the case of the Galileo I 
flyby a velocity increase of 3.92 mm/s was reported. This was a decrease of 
—2 mm/s for the Cassini spacecraft. Notice that, according to this model, 
the result would be similar for the most recent Juno flyby. 
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cannot attribnte a general validity to that equation. Using the gravitomag- 
netic model discussed in this paper we obtain AV ~ 0.65 mm/s for f3 = 10“^ 
in the last Rosetta flyby, in the same range of te prediction of Eq. (1), but 
different from the zero result obtained from the data analysis. 

We have also noticed that the value of /5 that we need to £t the anomaly in 
the case of the Cassini and Galileo I flybys is larger than the one correspond¬ 
ing to the other three flybys (NEAR, Galileo II and Rosetta). Inspecting 
Table 1 we can associate this inconsistency in the (3 values with the larger 
eccentricities in the Cassini and Galileo I orbits. On the other hand, we could 
elaborate an improved model in which f3 is not a constant but a function of 
r in order to obtain a better ht of the flybys. However, we are about to see 
that further speculation along this line is not required because the transver¬ 
sal gravitomagnetism phenomenological model cannot be made compatible 
with the flyby anomalies as well as the Gravity Probe B and the geodynamics 
satellites data. 


4 Observational Constraints 

A problem which arises in any putative model of the flyby anomaly is that 
we have already many constraints based upon careful measurements of the 
Earth’s gravitational held by using low orbit satellites. The extra perturba¬ 
tion detected in the hybys can be considered very large in the context of the 
present degree of accuracy achieved in these measurements. Consequently, 
we must discuss the contributions of any new interaction, postulated to ex¬ 
plain the hyby anomaly, in other experiments performed on satellites and 
check whether its contribution is sufhciently small to remain still undetected 
or not. There are two high-precision measurements of the Earth’s gravita¬ 
tional held (including gravitomagnetic effects): the geodynamics satellites 
and the Gravity Probe B experiment. 

The LAGEOS and LAGEOS II are spherical artihcial spacecraft entirely 
covered by retrorehectors which have been accurately tracked by laser rang¬ 
ing since their launch in 1976 and 1992, respectively. One of the objec¬ 
tives of these geodynamics satellites is to check the so-called Lense-Thirring 
effect, i. e., the precession of the longitude of the ascending node of the 
orbit as the consequence of the perturbation induced by the Earth’s gravit¬ 
omagnetic held as predicted by General Relativity. This ehect is very small 
(around 30 milliarcsec/year) in comparison with the much larger perturba- 


16 


tions contributed by the zonal harmonic coefficients of the Earth’s geoid 
[Ciufohni(1987)]. However, the classical nodal precession depends on the 
cosine of the inclination of the orbit in contrast with the relativistic Lense- 
Thirring effect which is independent of the orbital plane orientation. This 
gives a chance to deduce the tiny relativistic effect by comparing the data for 
several satellites. More recently, the similar LARES and GRACE missions 
added more data to the analysis. Despite these efforts, no conclusive results 
have been obtained yet [lorio et ah (2011)]. 

We should now calculate the contribution of the postulated extra gravit- 
omagnetic held in Eq. (5) to the precession of the longitude of the ascending 
node in these artihcial satellites. In terms of the true anomaly, z/, the equa¬ 
tions of the elliptic orbit can be written as follows [Valtonen and Karttunen(2005)]: 


r 


V 


cos n -h sin v s 
^1 -h e cos(z/ — iS) 



{ — (sin V -h esincu) h 


-h (cosz/-h ecoscn) s} , 


(17) 

(18) 
(19) 


where p = a(l — e^) is the semi-latus rectum and u is the argument of 
the pericenter, i. e., the angle between the radiovector of the periapsis and 
the direction of the ascending node. The system of reference in cartesian 
coordinates is given by h = i, s = cos I] + sin I k and w = cos J k — sin /j,/ 
being the inclination of the orbital plane. So h points towards the direction 
of the ascending node, w is the inclination vector and s forms a right-handed 
orthogonal system with them. 

The polar and azimuthal angles can also be expressed in terms of the true 
anomaly and the inclination of the orbit as follows: 


cos 6* = 

sin v sin I 

(20) 

cos 0 = 

cos u 

(21) 

\/l — sin^ / sin^ i/ 


The contribution of the perturbation forces to the instantaneous variation of 
the longitude of the ascending node, D, is given by [Burns (1976), Pollard(1966), 
Danby(1988)]: 


dVt 

dt 


r sin u 
77 sin J 
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where H = y/JI^ is the angular momentum of the orbiting spacecraft per 
unit mass and J\f is the component of the perturbing force perpendicular to 
the orbital plane. In our case, we have: 

Af = (3 Q— sin 9 cos 0 v x 0 ■ w , (23) 

r 

From Eqs. (17),(20) and (22) we finally have: 

R 

= dVL—. -|sin^ z/cos z/(sin J — cos/) 

a(l - e2) ^ ^ ^ (24) 

+ e cos uj sin^ z/ sin J — e sin u sin^ z/ cos z/ cos / } . 

Remarkably, every term in Eq. (24) cancels out when averaged over a full 
orbit, u G (0,27r). As the determination of the Lense-Thirring precession is 
plagued with difficulties, arising from the very large contributions of the 
Earth’s zonal harmonic coefficients [lorio et al.(2011), Ciufolini(1987)], it 
could be argued that a phenomenological law, as that in Eq. (5), could consis¬ 
tently fit the ffyby anomalies without impacting noticeably in the precession 
of the line of nodes of the satellites. However, if a sufficiently strong held of 
the form in Eq. (5) is assumed, the contribution to the precession of a gyro¬ 
scope in a low orbit exhibits large oscillations during the orbit which should 
have been found in the Gravity Probe B experiment [Everitt et ah (2011)]. 
The result shown in Fig. 5 implies a cancellation of the precession every half¬ 
orbit. Nevertheless, the average value of these oscillations is as large as 11.6 
arcsec for the (3 = 2x 10“^ parameter which we have shown to be compatible 
with the hyby anomalies. 

It is also interesting to provide an estimation of the value of the bounds 
on the value of (3 according to the Gravity Probe B experiment. The North- 
South orientation of the gyroscopes was measured with the average precision 
of 18.3 X 10“^ arcsec/yr in connection with the determination of the geode¬ 
tic effect which implied a secular reorientation of the gyroscopes in that 
direction [Everitt et ah (2011)]. Gonsequently, we can estimate a maximum 
non-secular average fluctuation of the North-South orientation of the gyro¬ 
scopes of this magnitude, which should correspond to < 3 x 10“®. Glearly 
insufficient to explain the flyby anomaly but, perhaps, within the range of 
other measurements in the Solar System. 


dVL 

dt 
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Figure 5: Contribution to the North-South orientation of a gyroscope in the 
Gravity Probe B orbit by a gravitomagnetic held of the form given in Eq. 
(5) with /3 = 2 X IQ-h 
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5 Conclusions and Remarks 


In this paper we have explored the idea of a strong gravitoniagnetic field 
as the sonrce of the flyby anomaly. Although this has been done in a pure 
phenomenological context is not at all excluded in principle. The condition 
of general covariance is too weak to restrict the form of the field equations of 
General Relativity, as many authors have recognized since the early days of 
the theory almost a century ago. Since them, many alternatives have been 
considered and dismissed [Goenner(2004)] and some with local and cosmo¬ 
logical implications, such as torsion models, are still a field of active research 
[Hehl(1976), Poplawski(2010)]. We have shown that a transversal component 
of a gravitomagnetic field would exhibit a maximum contribution to orbital 
perturbation effects in hyperbolic or highly eccentric elliptical orbits and it 
can be defined in such a way that the secular contribution to the Lense- 
Thirring precession is null. This is, precisely, what in the general features of 
the flyby anomaly is found. 

On the other hand, the detailed orbital evolution of the perturbations 
in gyroscope and Lense-Thirring precessions appear to be very large for a 
sufficiently strong field predicting a velocity variation in the millimeter range 
among the post-encounter and pre-encounter sections of the hyperbolic flyby 
orbit. Similar stringent constraints have been found in the earth-bound dark 
matter model [Adler(2009)]. The difficulties in settling interaction models 
verifying the conditions of being compatible with precise measurements of 
satellites in circular and elliptical orbits and also the large anomalies in the 
flyby data could point out towards a possible systematic error as the source 
of the discrepancies in the X-band Doppler residuals [Anderson et al.(2008), 
Acedo(2014)]. 

Anyway, many theoretical and computational work remains to be done 
in order to clarify all the conventional sources of perturbation on the flyby 
trajectories. In particular, a detailed numerical model of the atmospheric 
drag taking into account the density of the atmosphere and the geome¬ 
try of the spacecraft in a finite-element model is necessary to provide any 
reliable prediction of the dragging effects suffered by the spacecraft cross¬ 
ing the thermosphere during their closest approach to the Earth (this was 
specially noticeable in the case of the Galileo II flyby). The influence of 
albedo and infrared radiation will also be possible to analyze in detailed 
models of the spacecraft geometry as those build for the Pioneer 10 and 11 
[Rievers and Lannnerzahl(2011), Turyshev et al.(2012)]. 
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The solution to this riddle could also be found through the forthcoming 
Space-Time Explorer and Quantum Equivalence Principle Space Test (STE- 
QUEST) mission in which a spacecraft would be launched to a an elliptic 
orbit with very large eccentricity e = 0.773 specially suited to simulate an 
Earth flyby every period of 16 hours [Paramos and Hechenblaikner(2012)]. 
Another forthcoming mission: Gravity Recovery and Climate Experiment- 
Follow On (GRACE-FO) will also provide even more accurate data on the 
Earth’s gravitational held that could be instrumental in the next few year to 
elucidate this issue [Gruber (2010)]. 
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